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ABSTRACT. We show the existence of smooth stationary solutions for the inelastic Boltzmann 
equation under the thermalization induced by a host-medium with a fixed distribution. This is 
achieved by controlling the L p -norms, the moments and the regularity of the solutions for the 
Cauchy problem together with arguments related to a dynamical proof for the existence of stationary 
states. 



1. Introduction 

The dynamics of rapid granular flows is commonly modelled by a suitable modification of the 
Boltzmann equation for inelastic hard-spheres interacting through binary collisions |[T8l 1391 . As 
well-known, in absence of energy supply, inelastic hard spheres are cooling down and the energy 
continuously decreases in time. In particular, the Boltzmann collision operator for inelastic hard 
spheres does not exhibit any non trivial steady state. This is no more the case if the spheres are 
forced to interact with an external agency (thermostat) and, in such a case, the energy supply may 
lead to a non-equilibrium steady state. For such driven system (in a space homogeneous setting), 
the time-evolution of the one-particle distribution function f(v,t),v £ R 3 , t > satisfies the 
following 

d t f = rQ(f,f) + g(f), (1.1) 

where r ^ is a given constant, Q(f, /) is the inelastic Boltzmann collision operator, expressing 
the effect of binary collisions of particles, while Q{f) models the forcing term. 

There exist in the literature several physical possible choices for the forcing term Q in order 
to avoid the cooling of the granular gas: stochastic heating, particles heating or scaled variables 
to study the cooling of granular systems and even a nonlinear forcing term given by the quadratic 
elastic Boltzmann operator has been taken into account l25l . These options have been studied first 
in the case of inelastic Maxwell models Il2l[l9l[l3l[l4l[38l|6l|7l|2l][T6l. The most natural one is 
the pure diffusion thermal bath for which 

Q{j) = iiM (1.2) 

where fj, > is a constant, studied in ll26l[34l for hard-spheres. Such a forcing term corresponds 
to the physical situation in which granular beads receive random kicking in their velocity, like 
air-levitated disks [10]. Another example is the thermal bath with linear friction 



g(f)= li Af + \6xv(vf), 
i 



(1.3) 
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where A and fi are positive constants. We also have to mention the fundamental example of 
anti-drift forcing term which is related to the existence of self-similar solution to the inelastic 
Boltzmann equation: 

g(f) = -«div(«/), k>o. (i.4) 

This problem has been treated in |[3T1 |32l |33l for hard-spheres. For all the forcing terms given by 
(11.21) . (11.31 ). (11.41) it is possible to prove the existence of a non-trivial stationary state F ^ such 
that 

tQ{F,F)+G(F) =0. 

Moreover, such a stationary state can be chosen to be smooth, i.e. F £ C°°(IR 3 ). Finally, even if 
the uniqueness (in suitable class of functions) of such a stationary state is an open problem, it can 
be shown for all these models that, in the weakly elastic regime in which the restitution coefficient 
is close to unity, the stationary state is unique. For an exhaustive survey of the "state of art" on the 
mathematical results for the evolution of granular media see |39l . 

We are concerned here with a similar question when the forcing term Q is given by a linear 
scattering operator. This corresponds to a situation in which the system of inelastic hard spheres 
is immersed into a so called particles thermal bath, i. e. Q is given by a linear Boltzmann collision 
operator of the form: 

G(f) = B[f,F l ] 

where Fi stands for the distribution function of the host fluid and £>[•, •] is a given collision operator 
for (elastic or inelastic) hard-spheres. The precise definition of Q is given in Subsection 2.1. 

This kinetic model has already been tackled for instance in 13 in order to derive closed 
macroscopic equations for granular powders in a host medium. Let us also mention the work (51 
that investigates the case of a particles thermal bath made of elastic hard-spheres at thermodynam- 
ical equilibrium (i.e. Fi is a suitable Maxwellian). The deviations of the steady state (which is 
there assumed to exist) from the Gaussian state are analyzed numerically. For inelastic Maxwellian 
molecules, the existence of a steady state for a particles thermal bath has been obtained in |[2T1 . 
To our knowledge, the existence of a stationary solution of (11.11) for particles bath heating and 
inelastic hard- spheres is an open problem and it is the main aim of this paper. 

Our strategy, inspired by several works in the kinetic theory of granular gases ll26l l32l or 
for coagulation-fragmentation problems HEU, is based on a dynamic proof of the existence of 
stationary states, see I2T1 Lemma 7.3] for a review. The exact "fixed point theorem" used here is 
reported in Subsection 2.2. The identification of a suitable Banach space and of a convex subset 
that remains invariant during the evolution, is achieved by controlling moments and L p -norms 
of the solutions. In Section 3, we present regularity properties of the gain part of both collision 
operators Q and Q in (11.11) . Then, in Section 4 we get at first uniform bounds for the moments 
and the Lebesgue norms; in addition, we prove the strong continuity of the semi-group associated 
to (11.11) . and the existence and uniqueness of a solution to the Cauchy problem. All this material 
allows to obtain, in Section 5, existence of non-trivial stationary states. Finally, Section 6 contains 
the study of regularity of stationary solutions. Many technical estimates involving the quadratic 
dissipative operator Q(f, f) are based on results presented in l' T71[32l[33l[37l and in the references 
therein, but their extension to the linear inelastic operator G(f) is not trivial at all for the following 
reasons. First, since Q is not quadratic, it induces a lack of symmetry particularly relevant in the 
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study of propagation of LP norms. Second, since the microscopic collision mechanism is affected 
by the mass ratio of the two involved media (thermal bath and granular material), Povzner-like 
estimates for Q are not straightforward consequences of previous results from |[26l . Let us finally 
mention that our analysis also applies to linear scattering model which corresponds to the case 
r = 0. For such a linear Boltzmann operator, we obtain the existence of an equilibrium solution, 
generalizing the results of (30j [28] [38j to non-necessarily Maxwellian host distribution. 



2. Preliminaries 

Let us introduce the notations we shall use in the sequel. Throughout the paper we shall use 
the notation (•) = yl + | • | 2 - We denote, for any rj £ R, the Banach space 



measurable ; 



\f{v)\ (v^dv <+oo\. 



More generally we define the weighted Lebesgue space Lj)(R 3 ) (p £ [1, +oo), rj £ R) by the 
norm 

\f( V )\p(vrdv\ P . 



The weighted Sobolev space W v 
norm 



k > p 'T® 3 ) (p e [1, +oo), r] £ R and k £ N) is defined by the 

i/p 



w, 



k.p 



|s|<fc 



where d s D denotes the partial derivative associated with the multi-index s £ N . In the particular 

k 2 

by using the Fourier transform. 



case p = 2 we denote = Moreover this definition can be extended to for any s ^ 



2.1. The kinetic model. We assume the granular particles to be perfectly smooth hard spheres of 
mass m = 1 performing inelastic collisions. Recall that, as usual, the inelasticity of the collision 
mechanism is characterized by a single parameter, namely the coefficient of normal restitution 
< e < 1. To define the collision operator we write 

Q(fJ) = Q + (f,f)-Q-(fJ), 

where the "loss" term Q~(f, /) is 

Q-(f,f) = f(f*\v\), 

and the "gain" term Q + (/, /) is given by 

1 



(2.1) 
(2.2) 



Q + (fJ) 



Aire 2 



where the pre-collisional velocities read as 



>V = V + Ye ^ a 



\v — w\f('v)f('w) dadw, 

c 



w = w 



2e 



(\qW - <?), 



(2.3) 
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with ( = i^p. Notice that we always have ~ < £ < 1. Its weak formulation will be the main tool 
in the rest and it reads as 

/ Q + (f,f)(v)tp(v)dv = ±- ! [ f(v)f(w)\q\ [ ^(v')dadwdv, (2.4) 

,/R 3 4?r JR3 J R 3 J§2 

where q = v — w is the relative velocity of two particles about to collide, and v' is the velocity 
after the collision. The collision transformation that puts v and w into correspondence with the 
post-collisional velocities v' and w' can be expressed as follows: 

v' = v+-(\q\a-q), w' = w - | (|g|cr - g). (2.5) 

Combining (12.21 ) and (12.41 ) and using the symmetry that allows us to exchange v with w in the 
integrals we obtain the following symmetrized weak form 

/ Q{f,f)(y)i>{v)dv = \ I [ f(v)f(w)\q\A c ty)(v,w)dwdv, (2.6) 

./R 3 I J R 3 J R 3 

where 

Ac[ip](v,w) = — I N>{v') + Mw')-'4>{v)-'4>(w))d(j. (2.7) 

The inelastic Boltzmann operator Q(f, f) satisfies the basic conservation laws of mass and 
momentum, obtained by taking ifi = l,v in the weak formulation (12.61 ). since -4^[1] = -A^[v] = 0. 
On the other hand, in the modelling of dissipative kinetic equations, conservation of energy does 
not hold. In fact, we obtain A$ [\v\ 2 ] = — ^-r-\v — w\ 2 from which we deduce 

/ Q(f,f)(v)\v\ 2 dv = -^—^[ [ \v-w\ 3 f(v)f(w)dvdw, (2.8) 

JR 3 o JR 3 JR 3 

where we observe the dissipation of kinetic energy. In the absence of any other source of energy, 
the system cools down as t — ► oo following Haff 's law as proved in P2l . 

As already said in Introduction, the forcing term Q arising in the kinetic equation (11.11 ) is chosen 
to be a linear scattering operator, corresponding to the so called particles bath heating, 

G{f) :=£(/) = / / | 9 -n| [e- 2 /K)FiK)"/WFiH] dwdn (2.9) 
2tt\ J R 3 J § 2 

where A is the mean free path, q = v — w is the relative velocity, and w* are the pre-collisional 
velocities which result, respectively, in v and w after collision. The collision mechanism related 
to the linear scattering operator is characterized by 

(v — w) ■ n = —e{v± — w+) ■ n, (2.10) 

where n G S 2 is the unit vector in the direction of impact and < e < 1 is the constant restitution 
coefficient (possibly different from e). Here, we will consider a similar separation of the operator 
into gain and loss terms, C(f) = C + (f) — £~(f), with obvious definitions. Here the host fluid is 
made of hard-spheres of mass mi (possibly different from the traced particles mass m = 1) and 
the distribution function Fi of the host fluid fulfils the following: 
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Assumption 2.1. Fi is a nonnegative normalized distribution function with bulk velocity vl\ E 
and temperature Q\ > 0. Moreover, F\ is smooth in the following sense, 

Fi E Hg(R 3 ), Vs,5^0 



and of finite entropy / F\{v ) log F\(v) dv < oo. 



Remark 2.2. Notice that, since F\ E L3, itf of finite entropy, it is well-known [2, Lemma 4] that 
there exists some x > smc/j that 

u(v) := — !— I I \(v - w) ■ n|Fi(w)ckt;dn ^ xV 1 + M 2 Vu £ I 3 . (2.11) 

A particular choice of the distribution function Fx, corresponding to a host fluid at thermodynam- 
ical equilibrium, is the following Maxwellian distribution 

Notice however that our approach remains valid for more general distribution function. 

For particles of mass m = 1 colliding inelastically with particles of mass m\, the restitu- 
tion coefficient being constant, the expressions of the pre-collisional velocities (u*, w*) are given 
by dSHD 

1-/3 . . . 1-/3 . 

= u - 2a 1 _ 2 p W ■ n ) n > = w + 2(1 - a) - _ ^ (q • n) n, 

where a is the mass ratio and f3 denotes the inelasticity parameter 

mi . 1 — e 

a = — 1- E 0, 1 , = — - E 0, 1/2). 

1 + mi 2 

The post-collisional velocities are given by 

v* = v - 2a(l - (3) (q ■ n) n, = w + 2(1 - a)(l - /3) (g • n) n. (2.13) 

This linear operator can also be represented in a form closer to ( I2.6I ). By making use of the 
following identity lfTTll22l . 

(q ■ n) + (p(n(q ■ n)) dn = ^ J <p ^ \ da 

for any function ip, with q = q/\q\, we can rewrite the operator as 

= A / / l«l [e" 2 /(«*)FiK) - /(«)Fi(«j)] du;da (2.14) 

with 

1-/3 , , , s - , „ ^ 1-P , , , x 

«* = f - " x _ 2/? (9 - , w ir = w + {l- a ) l _ 2f3 (<? - \lW) ■ 

For such a description, the post-collisional velocities are 

v* = v - a(l - /3) (q - \q\a) , w* = w + (1 - a)(l - /3) (q - \q\a) . (2.15) 
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We consider Eq. dl.lb in the weak form: for any regular tp = t/j(v), one has 

/ f{v,t)^{v)&v= T - 
dt l M 3 2 



1 f(v,t)ip(v)dv = T - I I f(v,t)f(w,t) \q\A{[ip](y,w)dwdv 



+ T / / \q\f(v,t)F 1 (w)J e mv,w)dvdw (2.16) 



where 



JeMv, w) = ±-Jjq. n| - rP(v)) dn = -L / (, •( r ! ) - < ■( r) ) do: 



2.2. Proof of stationary states: basic tools and strategy. As stated in the Introduction, the final 
purpose of this paper is to prove the existence of a non-trivial regular stationary solution F ^ 
to (fTTTT ). Namely, we look for F G L 1 , F ^ such that 

rQ(F,F)+C?(F) = 0. (2.17) 

Remark 2.3. Notice that such a problem is trivial in the elastic case e = 1 and whenever F\ is the 
Maxwellian distribution (12- 12b - Indeed, in such a case, the Maxwellian equilibrium distribution 
M? of C provided by 11301 128[l38l is a stationary solution to di.il ) since QiM^M 9 ) = (elastic 
Botlzmann equation) and C(M^) = 0. 

The main ingredients are to show the existence of fixed points for the flow map at any time, 
and thus a continuity in time argument of the semi-group that allows to identify this one-parameter 
family of fixed points as a stationary point of the flow. Contraction estimates were used in (6l |2l] 
while in the hard-sphere case the Tykhonov Fixed Point Theorem was the tool needed 11261 124ll32l 

a. 

The exact result that will be used can be summarized as: 

Lemma 2.4 (Dynamic proof of stationary states). Let Y be a Banach space and (St)t^o be a 

continuous semi-group on Y such that 

i) there exists Z a nonempty convex and weakly (sequentially) compact subset of Y which is 
invariant under the action of St (that is StZ £ Z for any z £ Z and t ^ 0); 

ii) St is weakly (sequentially) continuous on Z for any t > 0. 

Then there exists zq G Z which is stationary under the action of St (that is StZo = Zofor any 

The strategy is therefore to identify a Banach space Y and a convex subset Z C Y in order to 
apply the above result. To do so, one shall prove that 

• for any f £ Y, there is a solution / € C ([0, oo), Y) to Eq. (fTTTT ) with f(t = 0) = / ; 

• the solution / is unique in Y and if /o G Z then f(t) G Z for any t ^ 0; 

• the set Z is (weakly sequentially) compactly embedded into Y; 

• solutions to (11.11) have to be (weakly sequentially) stable, i.e., for any sequence (f n ) n C 
C ([0, oo), Y) of solutions to (11.11 ) with f n (t) G Z for any t ^ 0, then, there is a subse- 
quence (fn k )k which converges weakly to some / G C ([0, oo), Y) such that / is a solution 

to CD- 
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If all the above points are satisfied by the evolution problem d 1 - lb - then one can apply Lemma 
to the semi-group (St)t^o which to any /o G Y associates the unique solution f(t) = Stfo 
to (ll.lb . Moreover, the regularity properties of the gain part of the operators [32] shall provide us 
the needed regularity to show the existence of smooth stationary states. 

3. Regularity of gain operators 



We recall the following result, taken from [32, Theorem 2.5, Proposition 2.6] and based on 
|fT7ll2"9l , on the regularity properties of the gain part operator Q + {g, /) that we state here only for 
hard-spheres interactions in space dimension N = 3. 

Proposition 3.1 (Regularity of the gain term Q + ). For all s, r\ > 0, we have 
\\Q + (g,f)\\ H s + i^C(s,r,,e) 



where the constant C(s, rj, e) > only depends on the restitution coefficient e G (0, 1], s and rj. 
Moreover, for any p G [1, oo) and 5 > 0, there exist 6 G (0, 1) and a constant C$ > 0, only 
depending on p, e and 5, such that 



[ Q + (f, f) f- 1 dv < C S \\f\\ p £~ 9) + 5 

JR 3 



P 

L i/p' 



On the other hand, the linear operator C(f) is quite similar to the quadratic Boltzmann operator 
associated to hard-spheres interactions and constant restitution coefficient e by fixing one of the 
distributions. In fact, it is possible to obtain the following similar result: 



Proposition 3.2 (Regularity of the gain term C + ). For all s,rj > 0, we have 



\jC+(f)\\ HS+1 <C( S ,T7,e) 



F ill^ +2 ll/ll^ +2 + l|Fi|| L i +2 H/|| L i +2 



(3.1) 



where the constant C(s, rj,e) > only depends on the restitution coefficient e G (0, 1], s and r]. 
Moreover, for any p G ( 1 , oo ) and 5 > 0, there exist q < p and a constant K$ > 0, only depending 
on p, e and 5, such that 



[ £+(/) f- 1 dv ^ K S IIPiIU, H/lli; 1 H/lUa + 5 (\\F 



l\\Ll U\\' L p 



+ l|Fi|U H/ILi H/ll^ 1 ). (3.2) 

Proof. The proof of these two estimates relies on the same steps as in Sections 2.2, 2.3 and 2.4 
of 11321 , see also ||37ll . We need just to have the same basic estimates as in their case. We start with 
the proof of (13.11) . An expression of the Fourier transform of C + can be obtained as: 

T [£+(/)] (0 := [ expHe ■ v)C+(f)(v) dv = -L / G(£ + ,£_)da 
Jr3 47rA J$2 

with G(v, w) = \v — w\f(v)Fi(w), G its Fourier transform with respect to (v, w) and 

e+ = (1 - a(l - 0))£ + a(l - 0)\Z\v, 6- = «(1 - P)£ - a(l - /3)|£|(7. 
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With this expression at hand, it is immediate to generalize to C + the regularity result in 
Theorem 2.5, Proposition 2.6] giving d3.ll ). 



Now, let us prove the second result. We first notice that, as in O, the gain operator C + admits 
an integral representation. Actually, even if it is assumed in [3] that Fi is given by the Maxwellian 
distribution (12.121) . a careful reading of the calculations of yields 



C- + f{v)=l f(w)k(v,w)dw, 



(3.3) 



where 

1 - 27 

2e 2 i 2 \v -w\ Jv 2 .(w-v)=o~ 1 V ' 2 7 



k(v,w) = n 2 , / F 1 (v + V2 + —z r L (w-v))dV 2 



with 7 = fty^ and 7 = (1 — a) 1 1 _ 2 / ^ - Arguing as in 11321 . we define the operator T related to 
the Radon transform: 

T : g€ L^IR 3 , dv) i-> T g{v) = — f g{fj,v + z) dz 

\V\ JzLv 

where fj, = 1 — ~ 2 ^ . For any h G M. 3 , let denote the translation operator T^fiv) = f(v — h), 
for any v G M 3 . Then, for any g G L 1 (M 3 , dv), one sees that 

(t w o T)(g)(v) = - — - — 1 / g(fi(v - w) + z) dz 



\V ~ W\ J z ±( v - W ) 

' g(v — w + z-\ — —(w — v) I dz, W, w G 



\V ~ W\ J z _L(v- w ) V 2 7 

Choosing g = T- w Fi leads to the following expression of the kernel k(v, w): 

1 



k(v, w) 



2e 2 7 2 



t w o T o r_ 



(Fi)(u), w,wG 



This previous computation is at the heart of the arguments of IT321 Theorem 2.2], from which 
one gets a version of Lions' Theorem ETl for a suitable regularized cut-off kernel with collision 
frequency of the form B m , n (\q\,q ■ a) = &s n (W\) &s m (<? ■ o), with $>s n smooth and with compact 
support [jj, n] , and bs m smooth and supported in [— 1 + — , 1 — — ] . More precisely, defining the 
smoothed-out operator in angular and radial variables Ct as in [32, Section 2.4]: 

C S mn (f) = T^[ [ BrnMq.^f^F^dwda (3.4) 

then, for any 77 G M + and any p > 1, there is C(p,r],m,n) > depending only on p, rj and 
(m, n), such that 

H4™,„(/)lk < Cfe^m^)!^!!!^ ll/ILi |n| (3-5) 
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for some q < p given by q = if p G (1, 6] while q = | if p G [6, +oo) (see IT321 Corol- 
lary 2.4]). In particular, Holder's inequality leads to 

p-i 

c^Wp-UvK ( f Fdv) P II^J^IlL^CK^ii/ii^iiFiii^ii/n^ 1 

\JR 3 J 

for some explicit constant C(m,n) > 0. 

Similarly, one can define the remainder part of C + which splits as 



r+ _ r+ _. r+ — £+ + /*+ + C + 

Jm,n Jim , n 11,0 m,n <J Am , n i\,£\,m , n 



with 



where $ fin (|g|) = |?| - $s„(M) and b Rm (q ■ a) = I - b Sm {q ■ a), q G R 3 ,a G S 2 . Holder's 
inequality provides 

JR3 i/p -i/P 

with p' such that - + A = 1, hence we have to estimate Lf, norms of Ct a , £~t Q , £pr> 

r P p ' ' OHm.n tibm.n titi m ,n 

for = — 1/p'. 

One can easily use ll32l Theorem 2.1] to prove that, for any 77 G R, 



for some explicit constant e(m) that, since the angular part of the collision kernel is such that 
linim^oo H&fymlUirs 2 ) = 0' converges to as m goes to infinity. 

It remains to estimate the norm of C^ s (/). We follow now the lines of l36l Chapter 9, 
p. 395] (which differs slightly from |[32l Proposition 2.6] and is more adapted to the linear case). 
Precisely, we split / as / = f r + f r c = /(w)X{|d<r} + f( v )X{\v\>r} f° r some r > 0. Then, as in 
IT361 p. 395], there is some positive constant C > such that 

II^C/OII^^^IIPill^JI/ll^ 

while 



ns m Jfrc)L^C— ||/|U HFilU 



|2+r,| + M "^i + n 

with A > 0. 
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Gathering all the above estimates we get, for 77 = —1/p', 

[ c^wr-idv^cWfWZ 1 (-HFiIM/IIl* + — 11/MiFiiu ) 

m ' n u i/ p \n 2 i/p r 2 i/p / 

IP I 111?. II _ II /II . II flip- 1 



+ e(m)(\\F 1 \\ Ll \\f\\ P LP +\\F 1 \\ LP ||/|| L x ||/||^ 
V i/p > p i/p, 

^(C-+e(m))\\F 1 \\ Ll \\f\f LP +(c— + e(m))\\F 1 \\ LP \\f\\ Ll . 
V n / 2 u i/p \ r J i/ P 2 ^ 1/p 

The proof follows then by choosing first m large enough then r large enough and subsequently n 
big enough. 

□ 

4. Regularity estimates for the Cauchy problem 

4. 1. Evolution of mean velocity and temperature. Let f(v, t) be a nonnegative solution to dl.lb . 
Define the mass density, the bulk velocity 

g(t)= I f(v,t)dv, u (t) = -^/ u/(t;,t)d« 

JR3 £>(t] ,/ R 3 

and the temperature 



Note that Eq. (12.161) for ip = 1 leads to the mass conservation identity £>(i) = i.e. 

g(t) = q(0) := 1. 

Now, Eq. (12.161 ) for ip(v) = v yields 

u(t) = - — ~^ / / |u-«;|(u -«;)/(«, tjPiCwJdudu;, Vt > 

which illustrates the fact that the bulk velocity is not conserved. To estimate the second order 
moment of /, let us introduce the auxiliary function: 



F(t)= / \v -w\ 2 f(v,t)F 1 (w)dvdw. 
Jr 3 Jm 3 

Notice that 

F(t)= [ \v-u 1 \ 2 f(v,t)dv + — Gi = 3e(t) + |u(£)- Ul | 2 + — 0i. (4.1) 
J R 3 mi mi 

In particular, to obtain uniform in time bounds of the mean velocity and the temperature, it is 
enough to provide uniform in time estimates of F(t). With the special choice ip(v) = \v — ui| 2 
one has 

-4<M(«, w) = C{1 ~J M {{fr-q- \q\) da = -((1 - C)k| 2 = "^M 2 
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while 

JeMv, w) = 2a 2 (l - (3) 2 \q\ 2 - 2a(l - (3)(q, v - m) 

= -2/e(l - K)\q\ 2 - 2n{q,w - m), i/,i»El 3 

q; 

with k = a(l — /3) = — (1 + e) £ (0, 1) and (•, •) denoting the scalar product. It is easy to see that 
F{t) = - {l ~f )T I I f(v,t)f(w,t)\q\ 3 dvdw 



2«(1 - k) 
A 




|u — u>| 3 /(f , i)Fi(io) dv dw 
^ f f \q\(q,VL 1 -w)f(v,t)F l (w)dvdw. (4.2) 

^ ,/R 3 ,/r 3 



/ /M) dv 

JR3 



Now, since / f(v, t)dv = 1 for any t ^ 0, Jensen's inequality yields 



f(w,t)\q\ 3 dw ^ 



v — / wf(w, t) dw 



3 



|« - u(i)| ; 



and consequently 

./'(/'. /)./'(»'. rlr chr ^ / |v-u(*)| 3 /M)d*; 

\ 3 / 2 / \ 3/2 



where we used again Jensen's inequality. In the same way, 



Qf Ju-u(t)| a /(«,*)du) = (se(t)) 



\ 3/2 

\v — w\ z f(v, t)Fi(w) dv dv 

Finally, the third integral in (14.21 ) is estimated as 



/ / \q\ 3 f(v,t)F 1 ( y w)dvdw > ( [ [ \v-w\ 2 f(v,t)F 1 ( y w)dvdvj] / =F(t) 3/2 . 

Jr 3 JR3 \Jr3 J R 3 J 



/ / \q\(q, ui — w)f(v, t)Fi(w) dv dw ^ / / \q\ 2 \ui — w\f(v,t)Fi(w) dv dw 

^2/ \v -u 1 \ 2 f(v,t)dv \w - Ui|Fi(«;) dw + 2 / |io - ui| 3 Fi(u;) dw 

Jr 3 Jr 3 jr 3 



^ C F(t) 

where 

Cn = 2 max < 



/ | to — ui| 3 Fi(u>) dw 

/ |u>-ui|Fi(w)dw, ^ 3 
Jr 3 / 

JR 



to — ui| 2 Fi(u>) dw 
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In conclusion, we obtain 

Ht) < -^^(3e(t)f 2 - 2K{1 - K) F(t)^+ 2 ^F(t) < - 7 iF(t) 3 / 2 + 72 F(t) (4.3) 

2«(1 - «) , 2C k . . , „ , . . . , , 

where 71 = > and 72 = — - — > 0. A simple use of the maximum principle shows 

A A 

that 

F(i)^maxj^^ ,^(0)|, Vi > 0. 

Because of (14.11 ), this leads to explicit upper bounds of the temperature 0(f) and the velocity 
|u(t) — Ux|, namely 

sup ^39(t) + |u(t) - ui| 2 ^ ^maxH^ ,F(0)J<oo. (4.4) 

4.2. Propagation of moments. To extend the previous basic estimates, in the the spirit of lfl5l . 
we deduce from Povzner-like estimates some useful inequalities on the moments 



Y P (t)= / f(v,t)\v\ 2r dv, £>0, r^l 
in 3 

where f(t) is a solution to (11.11 ) with unit mass. One sees from (11.11 ) that 

—Y r (t) = rQ r (t)+L r (t), 
at 

where 

Qr{t)= I Q(f,f)(v,t)\v\ 2r dv, L r (t)= [ £(f)(v,t)\v\ 2r dv. 

The calculations provided in 1126*1 [T5l allow to estimate, in an almost optimal way, the quantity Q r . 
One has to do the same for L r (t) given by 

L r{t) = \ I I f(v,t)F 1 (w)\v- W \J e [\-\ 2r ](v,w)dvdw. 

A Jr 3 Jr 3 

To do so, let us derive Povzner-like estimates for C in the spirit of lf26l . The application of 
the result of [26 ] is not straightforward since, obviously, C is not quadratic and because of the 
influence of the mass ratio a = in the collision mechanism. Here, we will write the mass 

m+mi ' 

of particles m even if taken as unity for the sake of the reader. To be precise, we are looking for 
estimates of 

Je[\-\ 2r ](v,w) = ±- [ \q-n\(\v*\ 2r -\v\ 2r ) dn, r > 1. 

To do so, it shall be convenient to write 

Je[\ • \ 2r ](v,w) = — *— f \q-n\\ty (m\v*\ 2 ) - ^ (m\v\ 2 )} dn (4.5) 
27rm r Jg2 k v ' 

where ty(x) = x r , r ^ 1. We adopt the strategy used in lf26l and write 

^ (m\v*\ 2 ) - ^ (m|t;| 2 ) = q e (^)(v,w) + ^(mi|u;| 2 ) - ^(mi\w*\ 2 ) (4.6) 
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where 

q e (^)(v,w) = * [m\v*\ 2 ) + ^ (mi\w*\ 2 ) - ^ (m\v\ 2 ) - ^ (mi\w\ 2 ) . 

Now, 

q e (*)(«, W) = p e (* )(«, ™) - tle(*)(«, W) 

with 

{p e (^')(v, to) = * (m|t;| 2 + mi|to| 2 ) — * (m|-u| 2 ) — ^ (mi|w;| 2 ) 
n e (*)(«,iy) = * (m|u| 2 + mi|w| 2 ) - ^ (m|t?*| 2 ) - ^ (mi|w*| 2 ) . 

Applying ||26j Lemma 3.1] to the function ^ with x = m\v\ 2 and y = mi\w\ 2 , we see that there 
exists A > such that 

p e (*)(t», to) ^ A^m|u| 2 ^' (mi\w\ 2 ) + mi|u;| 2l I' / (m|w| 2 ) J (4.7) 

while, since \& is nondecreasing and m\v\ 2 + mi|u>| 2 ^ m|u*| 2 + mi|u)*| 2 , there exists b > 
such that 

n e (*)(u, to) ^ 6mmi|t;*| 2 \w*\ 2 $>" (m|-u*| 2 + mi\w*\ 2 ) . 

One can then write 

n e (^)(v,w) ^ b A(v*, to*) (m\v*\ 2 + mi\w*\ 2 ) 2 ^" (m|t;*| 2 + mi\w*\ 2 ) 

where 

* *\ m|tJ*| 2 mill/;*! 2 

A{v ,w ) = g • 

(m|t;*| 2 + mi|to*| 2 ) 

To estimate better the above term A(v*, to*), it will be convenient to parametrize the post-collisional 
velocities in the center of mass-relative velocity variables, which, with respect to the usual trans- 
formation (see e.g. E6l Eq. (3.10)]) depend on the masses m and mi. Namely, let us set 

. z + m\t\q\m z — ml\q\vo 

V = , w = 

m + mi m + m\ 

where z = mv + miw,q = v — w and w is a parameter vector on the sphere S 2 . The parameter t is 
positive and such that v* — w* = £\v — w\dj. In particular, one sees from the representation (12.131 ) 
that < I ^ 1 . In this representation, one has 

|f*| 2 = t ttt ( U| 2 + m\i 2 \q\ 2 + 2^mi|g||z| cos if 

(m + m\Y \ 

and 

|to*| 2 = t wt| U| 2 +m 2 £ 2 |g| 2 - 2£m|g||,dcos/z ), 

(m + mi) 2 \ 

where /z is the angle between z and tu. One has then 

m\v*\ 2 + mxltt;*! 2 = ( \z\ 2 + £ 2 mmi\q\ 2 ) . (4.8) 

m + mi 
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One can check that 

(m|i;*| 2 ) (mi|u;*| 2 



\z\ + 1 mmi\q\ 



+ 



(?7l + TOi) 4 

£(toi — m)\z\ \q\ + (\z\ 2 — £ 2 mmi\q\ 2 ^j 



z\ —I TOTOll^l 



cos 2 p 



COS [I 



4£ 2 mmi\z\ 2 \q\ 2 cos 2 p 



i.e. 



Therefore 



(m|t;*| 2 ) {mi\w*\ 2 ) ^ 



Til nil 



(TO + TOi) 4 - 



|z| 2 + £ 2 mmi\q\ 2 (1 — cos 2 



. , . . . to mi n 
A(v*,w*) ^ — ; , 9 sin 2 /i. 



(TO + TOi) 2 

We obtain then an estimate similar to the one obtained in l26l . Moreover, it is easy to see from (14.8 
that 

m|v*| 2 + mi|ui*| 2 ^ I 2 (m\v | 2 + to-i|u)| 2 ) 
and, arguing as in [26], there exists some constant n > such that 



n e [^ , ](v, w) ^ sin 2 (m\v | 2 + mi|u;| 2 ) 2 (m|u| 2 + mi|u;| 2 ) . 



(4.9) 



This allows to prove the following: 



Lemma 4.1 (Povzner-like estimates for C). Let ^f(x) = x r , r > 1. Then, there exist positive 
constants k r and A r such that 



\v-w\J e [\ ■ \ 2r ](v,w) < ,4 r (jw|H 2r + |v| 2r |w|J + ^p\v - w\\ir\-' - /,> ( r 
/or arcy u, u> G M 3 . 

Proof. Bearing in mind that J e \\ ■ \ 2r ](v,w) is provided by (14.51 ) and (14.6b . first of all, since 
^ (toi|ui*| 2 ) ^ 0, we note that 



2r 



2r+l + , ,2r+l 



^(toIv^I 2 ) - ^J/(m|t;| 2 ) ^ q e [ty](v,w) + * (toi|u;| 2 ) = q e [^](v, w) + 



r i „ , , 1 2r 



IP 



Then, integrating (14.71 ) and ( 14.91 ) with respect to the angle n G S 2 , one obtains, as in |[26l Lemma 
3.3.] and [26, Lemma 3.4], that there are A r and k r > such that, for any v, w G M 3 : 



m — w\ 



l — jf q e (*)(«,«;) | g • n| dn < A r (\v\\w\ 2r + |u | 2r |u»|^ - k r (\v\ 2r+l + \w\ 2r+1 ^j . 



27TTO. 

and this concludes the proof. 



□ 



The above Lemma (restoring m = 1) together with the known estimates for Q r (t) allow to 
formulate the following 

Proposition 4.2 (Propagation of moments). Let f(t) be a solution to (11.11 ) with unit mass. For 
any r ^ 1, let 

Y r (t)= [ f(v,t)\v\ 2r dv, t^O. 
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Then, there are positive constants A r , K r and C r that depend only on r, a, (3, r, A and the mo- 
ments o/Fi such that 

-r^r(t) ^ C r + A r Y r (t) - K r .Y r 1+1 / 2r (t), Vt ^ 0. 

As a consequence, if Y r (0) < oo, then sup t ^ Y r (t) < oo. 
Proof. Recall that -$f¥ r (t) = rQ r (t) + L r (t), where 

' 2r j„, T U\ — I rf f\f„, +M„.|2r . 



Qr(*)=/ Q(/,/)(^t)H" r d^ ; L r (t)=/ £(/)(M)|< r d<;. 

According to |[26l Lemma 3.4], there exist A r > and fc r > such that 

Q r (t) ^ A r Y 1/2 {t)Y r (t) - k r Y r+1/2 (t), t ^ 0. 



Now, from Lemma |4~T1 



XL r (t) sC ^ r M r Y 1/2 (t) + A r M 1/2 Y r (i) +ml \ 



|u — w\\w\ 2r f(v ,t)Fi(w) dw dv 
- k r Y r+ y 2 (t) - k r M r+ y 2 , 



where M s = dw, s ^ 1. One has 

/ / |t;-t(;||u;| 2r /(w,t)Fi(t(;)d'u;du^M r Y 1/ 2(t) + M r+1/ 2 
and, denoting c 1 / 2 '■= sup t ^ Y 1 / 2 (t) < oo, one has 

L r (t) < C r + -JLJJlY r (t) - -f Y r+1/2 (t) 

where C r = [c\j 2 A r M r + c\j 2 m\M r + m^M r+1 / 2 ) /A is a positive constant depending only on 
a,/3,X, r ^ 1 and the moments of Fi. Gathering all these estimates leads to 

^Y r (t) < C r + A r Y r (t) - K r Y r+1/2 (t) 

where A r = rA r Ci/ 2 + jA r M 1 / 2 > and K r = r&v + ^ > 0. Now, thanks to the mass 

conservation and Holder's inequality, one gets Y r+1 / 2 (i) ^ Yl +1 ^ 2r (t) which leads to the desired 
result. □ 

Remark 4.3. We see from the definition of the positive constants A r , C r and K r that the above 
Proposition still holds true whenever r = (i.e. for the linear problem). 
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4.3. Propagation of Lebesgue norms. Let us consider now an initial condition /q G L\ n L p for 

some 1 < p < oo. We compute the time derivative of the L p norm of the solution f(v, t) to dl.ll) : 



l_d 



JR 3 J R 3 J R 3 

+ / C + (f)fv- l &v- [ C-{f)r-Ui 

JR 3 JRZ 



Using the fact that / f p ~ 1 Q~~(f,f)dv ^ and £"(/)(«) = v{v)f{v) where the collision 
frequency v(v) is given by 



u ( v ) = — r / / \(v — w) ■ n\Fx(w) dw dn, 
2vrA J R 3 J § 2 



we obtain the estimate: 



~f f(v,t)dv^T [ Q + (f,f)f p - 1 dv+ f £+(f)f- l dv- [ v{v)f p {v,t)dv. 

P dt J R 3 J R 3 J R 3 J R 3 



Using the lower bound (12- lib , we get 
l_d 



■\\f{t)\\%^ r I Q + {fJ)f- 1 dv+ I _C + (f)p- 1 dv- X \\f\\ P LP 



Proposition 13.11 and the conservation of mass imply that, for any 5 > 0, there is 6 > and some 
Cs such that 

/ Q + (fj)f p -Hv,t)dvKc 5 \\fm p £~ e) +wmLi\\fmi* ■ 

Moreover, Proposition |3.2| implies that, for any 5 > 0, 

c+fafip-^v^awmwF + <h& fn/(t)ii^ +n/(t)ii L i||/(t)ii^ 1 ), 

\ i/p i/p/ 

for some constants Ci, C2 > that depend only on p, 5, 77, a, e and the norms of Fi in the spaces 
involved in (13.21) . Recall that there is some M2 such that 



sup||/(i)|U = 1 + sup / |u| 2 /(«,*)dtJ < M 2 < 00. 
t^o 2 t^o Jm.3 

Now, using Young's inequality, xy p ~ l ^ -x p + p —^-y p , for any 2, y ^ 0, we have 



/ C + (f,f)f p - 1 dv < cy/WH^ 1 + C 3 <5 f 
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for some constant C3 > 0. Collecting all the bounds above, we get the estimate 

^rC 5 \\f(t)\\ P S~ 9) + 5rM 2 \\f(t)\\ P LP + C x \\f(t)\\l^ 

y CL£ 1/p 



+ C 3 5 (||/(t)||^ + ^- x 



1/p 



<;TC & \\f{t)\\ P £~ 6) + (W + C7 3 ) - |) ||/(t)||^ 



X 



V 



+ C 1 \\f(t)\\ p LP 1 + C 3 5MV-- UJy , JULP , 

since || • \\ L v ^ || • \\lp. Choosing now 5 such that 5(tM2 + C3) < x/2, we get the existence of 

i/p 

positive constants C4, C5 and Cq such that 

"^ll/Wlli* < C 4 ||/(t)||S- e) + ^/(i)^ 1 - C 6 \\f(t)\\ p LP + C7 3 M|5. 
It is not difficult to get then that sup f ^ ||/(i)||z,p < 00. This can be summarized in the following 

Proposition 4.4 (Propagation of L p -norms). Let p G (1, 00) and /o G L\r\ L p with unit mass. 
Then, the solution f(t) to (11.11 ) satisfies the following uniform bounds 

aup(||/(t)|U + ||/(t)||ip) <oo. 

Remark 4.5. Notice that the fact that Fx is of finite entropy (see Assumption \2. i| ) /jos teen Msec? 
/jere above, via the lower bound ( 12.111 ), in order to control from below L p norms involving the 
loss operator C~. Whenever r > 0, if is possible then to replace such estimates involving £~ by 
others that involve Q~. Notice also that, whenever r = (i.e. in the linear case), only the above 
constant C4 vanishes and we still have sup t ^ < °°- 

As a corollary, we deduce as in [32, Section 3.4], see also EOl . the following non-concentration 
result: 

Proposition 4.6 (Uniform non-concentration). Let /q be given with unit mass. Assume that there 
exists some p G (1, 00) such that /o G L\f\ L p . Then, there exists some positive constant Vq such 
that 

H)< f \v-u(t)\ 2 f(v,t)dv < 1/i/q, Vi^O, 

where f(v, t) is the solution to (11.11 ) with /(0) = /o awcf u(i) = / vf(v, t) dv, t ^ 0. 

Proof. Let /(i) be the solution to (11.11 ) with /(0) = /o- From the above Proposition, there exists 
C p > such that sup t ^ ^ Cp> an d Holder's inequality implies that, for any r > 0, 



sup/ /(M)dw<C p — r 3 

t^O J{|»-u(i)|<r} V <3 
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Accordingly, there is some ro > such that 

f(v,t)dv^l, Vi^O. 

f>-u(t)|<r } 1 

Then, for any t ^ 0, recalling that / f(v, t) dv = 1 for any t ^ 0, 

f(v,t)\v-u(t)\ 2 dv> f f(v,t)\v-u(t)\ 2 dv^r 2 f f(v,t)dv 



>r 2 (l- / /(«,*) dt; >^ 

\ J{|^-u(t)|<r } / 2 



which concludes the proof. □ 

4.4. L 1 - stability. As in 11321 . in order to prove the strong continuity of the semi-group (St)t^o 
associated to (11.11 ). one has to provide an estimate of \\f(t) — g(t)\\ for two solutions f(t) and 
g(t) of (11.11 ) with initial conditions /(0),g , (0) in some subspace of L . This is the object of the 
following stability result, inspired by |[32l Proposition 3.2] and OTl Proposition 3.4]. 

Proposition 4.7 (L l -stability). Let /o, go be two nonnegative functions of L\ and let f(t),g(t) G 
C(M + ,L2) n L°° (M + , Lg) &e associated solutions to (ll.ll ). 77je«, « A > depending 
only on sup t ^ + such that 

Wfit) - g(t)\\ L i < ll/o - flttkj exp(At), Vt > 0. 

Proof. Let = /(i) — g(t). Then, /i satisfies the following equation: 

d t h(v,t)=T!.Q(f,f)-Q(g,g)) +£(h), h(0) = f -g . (4.10) 



As in (3T][32l, the proof consists in multiplying (14.101 ) by ip(v, t) = sga(h(v, t))(v) 2 and integrat- 
ing over M 3 . We get 



d 

di 

where 



\h{v,t)\(vydv = I(t)+L(t) 



I(t)=rf {Q(f,f)-Q(g,g))^(v,t)dv and L(i) = / C(h)(v,t)tf>(v,t)dv. 



To estimate the integral I(t) we resume the arguments of [31 , Proposition 3.4] that we shall need 
again later. According to ( 12.61 ) 



= I [ [ CfO, t) f(w, t) - 9(v, t)g(w, t)) \q\A c mt)](v, w) dw dv. 

2 7R3 J R 3 

The change of variables (v,w) i— > (w,v) implies that 



T 



/ W = 2 /- I „(f(v,t) ~ g(v,t)) (f(w,t) + g(w,t)) \q\A c mt)](v,w) dw dv. 
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Moreover, it is easily seen from the definition of ip that 

(f(v,t) - g(v,t))Acm)]M < ^ \fM - g(v,t)\ ({v') 2 + (w') 2 - {v) 2 + (w) 2 ) da 

^2\f(v,t)-g(v,t)\(w) 2 
where we used the fact that A([(-) 2 ](v, w) = — ^-j-\q\ 2 ^ 0. Therefore, 

Ht)<-r f [ \q\\f(v,t)-g(v,t)\(f(w,t)+g(w,t))(w) 2 dvdw 
Jr 3 Jr 3 

< r / \f(v, t) - g(v, t)\ (v) 2 dv [ (f(w, t) + g(w, t)) (w) 3 dw 

JR 3 JR 3 

i.e. 

I(t)^r\\f(t)+g(t)\\ Ll \\f(t)-g(t)\\ Lh , Vt > 0. (4.11) 

On the other hand, recalling that C(h)(v, t) = C + (h)(v, t) — i/(v)h(v, t), from formula ( 13.31 ) one 
has 

L(t)= I tp{v,t)dv / h(w,t)k(v,w)dw - / u(v)\h(v, t)\ (v) 2 dv 
Jr 3 Jr 3 Jr 3 

</ (v) 2 dv \h(w,t)\k(v,w)dw - u(v)\h(v, t)\ (v) 2 dv, 
Jr 3 Jr 3 Jr 3 

i.e. L(t) ^ / £(\h\)(v, t)(v) 2 dv. Now, since / C(\h\)(v, t) dv = for any h, one gets that 

Jr 3 Jr 3 

L(t) s$ f C(\h\)(v,t) \v\ 2 dv. 



Resuming the calculations performed in Section |4~T1 (see Eq. (14.21 )), one gets that 



C(\h\)(v,t)\v\ 2 dv < - 2/t(1 K) / / Iv-wflhiv^^F^dvdw 

* Jr 3 Jr 3 

+ ~T f f \q\(q,-w)\h(v,t)\Fi(w)dvdw 



X 



< y / / M 2 HK'f'--')iFt(«-)<lrd/r. 



This leads to 



L(i)<-^(2/ \v\ 2 \h(v,t)\dv [ \w\Fi(w)dw + 2 [ \h(v,t)\dv[ H 3 Fi(u>) dw 
X V ./r 3 in 3 Jr 3 Jr 3 

and, setting c+ = max { / \w\Fi(w)dw, f \w\ 3 Fi(w) dw 1 , we get 
L(t)^c+ / |%,t)|( ? ;) 2 dt; = C+ ||/(t)- 9 (t)|| Ll . 
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Gathering ( 14.1 11) together with the latter estimate and denoting then A = r sup t ^ \\f(t)+g(t)\\ L i + 
c+, we get the estimate 

■^ll/(*)-^)IU< A ll/(*)-^(*)ILi 5 

and the proof is achieved. □ 

4.5. Well-posedness of the Cauchy problem. We are in position to prove that the Boltzmann 
equation (11.11) admits a unique regular solution in the following sense: 

Theorem 4.8 (Existence and uniqueness of solution to the Cauchy problem). Take an initial da- 
tum /o G L\. Then, for allT > 0, there exists a unique solution f G C([0, T\; L\) n L°°(0, T; L\) 
to the Boltzmann equation f li.il ) smc/z f/ja? f(v,0) = fo(v). 

Proof Let T > be fixed. The uniqueness in C([0,T]; L\) n L°°(0, T; L|) trivially follows from 
Proposition 14.71 The proof of the existence is made in several steps, following the lines of OT1 
Section 3.3], see also 0511251 . 

Step 1. Let us first consider an initial datum /o G L\, and define the "truncated" collision 
operators 

Qn(f,f)(v)ip(v)dv = l f [ l M<n} \q\f(v)f(w)A C [i>](v,w)dwdv, 

2Jr3Jr3 (4.12) 

C n (f)(v)ip(v)dv = - / l{|g|< n } \q\f(v)F 1 (w)J e [tl)](v,w)dvdw 
A Jr 3 Jm 3 

for any regular test function tp. The operators Q n and C n are bounded in any L x q , and they are 
Lipschitz in L\ on any bounded subset of L\. Therefore, following (T), we can use the Banach 
fixed point theorem to get the existence of a solution ^ /„ G C([0, Tj; L\) nL°°(0, T; L|) to the 
Boltzmann equation 9t/ = rQ n (f, f) + C n {f). Thanks to the uniform propagation of moments 
in Proposition 14.21 there exists a constant Cr > (that does not depend on n) such that 

sup || f n \\ L i ^ Cr, VnGN. 



[0,31 



Step 2. Let us prove that the sequence (f n ) n is a Cauchy sequence inC([0, T]; L2)nL°°(0, T; L4). 
For any m ^ n, writing down the equation satisfied by f m — f n and multiplying it by ip(v , t) = 
sgn(/ m (v , t) — f n (v, t))(v} 2 as in the proof of Proposition 14.71 we get 

/ \f n (v,t) - f m (v,t)\(v) 2 dv = I m ,n{t) + Jm,n{t) 
JR 3 



dt 



where 



and 



(t) =T \Q m {fmJ m ) ~ Qn(fn,fn)\i>(v,t) dv 
Jm,n(t) = / \ C m (f m )(v,t) - C n (f n )(v,t) \lp(v,t)dv. 
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We begin by estimating I m ,n{t)- It is easy to see that I m ,n{t) = 4ra(*) + n(*) where 
J m,n(*) = 77 / / (fm(v,t) f m (w,t) - fn(v,t)f n (w,t)) B m (q)A^[ip(t)](v,w) dw dv, 

2 JR 3 J R 3 V / 



while 



Il,n(t) = I / / (B m (q) - B n (q)) f n (v,t)f n (w,t)A C {i>(t)](v,w) dw dv, 
1 Jr 3 Jr 3 

where B n {q) = Ml{|g|<n}- Arguing as in the proof of (14. Ill ), we get easily that 

4,n(*) < r\\f n (t) + / ro (t)|| L i||/„(t) - /m(*)|| £ i, Vt > 0. 

The estimate of 1^ n (t) is more involved. One observes first that 

B m (q) - B n (q) = |g|l{ n ^| g |<; m } < \q\ (l{|„|^ n /2} + l{H>n/2}) • 
As in the proof of Proposition 14. 7 1 one has 

^(t)](«,«0<T- / ((^) 2 + M 2 + ^) 2 + H 2 ) da^2((t;) 2 + H 2 ) 



4vr _ 

and, since \q\ ^ (v}(w), one gets 

4,n(*)^ / / /„(«,t)/«K*)|g|(l {\v\>n/2} + !{|«>|>n/2}) ((^) 2 + (w) 2 ) dw dv 
JR 3 JR 3 

/ fn(v,t)f n (w,t){v){w) ((vf + (wf) (l{/„)W2} + l{(w)>n/2}) dwdv. 
JR 3 JR 3 

It is not difficult to deduce then that 

4,n(*) ^ 4r (^ 3 ^'^^{j^ /»( U >*)< W > 3l {<«W2} d ^ 

Since sup[ T ] H/nC^Hi, 1 ^ Ct for any n G N, the latter integral is estimated as 



/n(M)(«) d l{>)>n/2} < / / n (u,t)(u) 4 l {(t , )>n/ 2}-^ ^ — ^ 

\v) n 



and we get 



VJr 3 ) n n 

Therefore, 

8r 2 r 

/m,n(t)<2rCr||/ n (t)-/m(t)Li + — I -, Vie[0,T], m^n. (4.13) 

We proceed in the same way with J m>n {t). First, we notice that J m , n {t) splits as J m>n (t) = 
4,n(*) + 4,nW with 

J 1 m,n( t ) = \l I B m (q)[f m (v,t)-f n (v,t)]F 1 (w)J e [iP(t)}(v,w)dvdw 

A ./IBS ./ib3 
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and 

= T / f i B m{q)-B n (q)]f n (v,t)F 1 (w)J e [ij(t)](v,w)dvdw. 

A JR 3 Jr 3 

Arguing as in the proof of Proposition 14.71 we get 

/ ^m(|/n-/m|)(«,t)<«> 2 d«<^ / / 1 ?| >| | (/„ " fm) (v, t) \F, (w) dv dw 

and there exists a positive constant c+ such that 

J^ n (t) < c+||/ n (i) - / m (t)|| L i, Vi £ [0,T]. 
Let us now estimate J^ n (t). As above, 

J mA t ) = Tl I \ ( l\ 1 {n^\q\^m}fn(v,t)Fi(w)J e [i;(t)}(v,w)dvdw 

and 

J e m)]M<^- f \q-n\((v*) 2 + (v) 2 ) dn = J e [(.) 2 ]+^-[ \q ■ n\(v) 2 dn. 

Calculations already performed lead then to 

J e [t/j(t)](v,w) ^-2 K (q,w) +2(v} 2 ^2((v)(w} 2 + (v) 2 ) , V«,»GR 3 . 

Finally, 

<»(*) < t / / k|i{| 9 |^}/n(^,t)FiH ((^)M 2 + («) 2 ) dwdv 

A JR 3 JR 3 

< 7 / / l{k|^}/n(^,t)FiH ((«) 3 M + (v) 2 (w) 3 ) dwdv. 
A Jr 3 Jr 3 

Now, arguing as we did for I^^ n {t), there exists some constant Ct that depends only on ||Fi|| L i 
and sup n sup [0 T] ||/ n (*)||xi such that 

— > VrrOn, ie[0,T]. 
n 

Gathering all these estimates, we obtain the existence of constants C\(T) and C2(T) that do not 
depend onm,ti such that 

4l I |/n(^,t)-/ m («,t)|(«) 2 d V ^C 1 (T)||/ n (t)-/ m (t)|| Ll + ^) ; vte [0,T],m>n. 
at Vina Mij 2 n 



This ensures that (/ n ) n is a Cauchy sequence in C([0,T]; L^). Denoting by / its limit, we 
obtain that / £ C([0, T]; Z^) n L°°(0, T; L|) is a solution to the Boltzmann equation (fTTTT) (with 
the actual collision operators Q and C). 

Step 3. When the initial datum /q £ Lg, we introduce the sequence of initial data foj := 
/ol|u|<j- Since foj £ L4, we have the existence of a solution fj £ C([0, T]; L\) n L°°(0, T; L4) 
to the Boltzmann equation associated to the initial datum foj. Moreover, there exists Ct such 
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that supp T ] ||/j||x,i ^ Cr- We establish again that {fj)j is a Cauchy sequence in C([0, T]; L\) by 
using the L 1 -stability in Proposition 14.71 □ 

5. Existence of non-trivial stationary state 

All the material of the previous sections allows us to state our main result: 

Theorem 5.1 {Existence of stationary solutions). For any distribution function ~F\{v) satisfying 
Assumption \2. 1 1 and any r ^ 0, there exists a nonnegative F £ L^H LP, p £ (1, do) with unit 
mass and positive temperature such that tQ(F, F) + C{F) = 0. 

Proof. As already announced, the existence of stationary solution to (11.11 ) relies on the application 
of Lemma 12.41 to the evolution semi-group (St)t^a governing (ll.ll ). Namely, for / £ L l , let 
f(t) = Stfo denote the unique solution to (11.11 ) with initial state /(0) = /o The continuity 
properties of the semi-group are proved by the study of the Cauchy problem, recalled in Section[4] 
Let us fix p E (1, oo). On the Banach space y = L\, thanks to the uniform bounds on the L\ and 
LP" norms, the nonempty convex subset 

Z = ^fey, J fdv = l and \\f\\ L i + ||/|| iPo < m\ 

is stable by the semi-group provided M is big enough. This set is weakly compact in y by 
Dunford-Pettis Theorem, and the continuity of St for all t ^ on Z follows from Proposition [477] 
Then, Lemma 12.41 shows that there exists a nonnegative stationary solution to (11.11 ) in L\ D IP° 
with unit mass. In fact, the uniform in time LP bounds also imply the boundedness of F in LP for 

ah>e(l,oo). □ 

As a corollary of Theorem 15.11 choosing r = allows us to prove the existence of a steady 
state to the linear inelastic scattering operator C when the distribution function of the background 
is not a Maxwellian, generalizing the result of Il30ll2"8ll38l . 

Corollary 5.2. Let Fi satisfy Assumption \2.1\ Then, the linear inelastic scattering operator C 
defined by (12.91 ) admits a unique nonnegative steady state F G L\v\ L p , p £ (1, oo), with unit 
mass and positive temperature. 

Proof. The existence of a nonnegative equilibrium solution F £ L\ is a direct application of 
Theorem [5J] with r = 0. A simple application of Krein-Rutman Theorem implies the uniqueness 
of the stationary solution F within the range of nonnegative distributions with unit mass. □ 

Remark 5.3 (H-Theorem and trend towards equilibrium). As in [28], it is possible to prove a 
linear version of the classical H-Theorem for the linear inelastic Boltzmann equation (11.11 ) with 
T = : 

dtf = £(/), /(t = 0) = / £L 1 . (5.1) 
Namely, for any convex C 1 function $ : R + — > R, let 

H*(f\F) = [ F(v)$ ($\) dv, f £ L\ 
J R3 \F(v)J 
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Arguing as in ff28i it is easy to prove that, if the initial state /o has unique mass and finite entropy 
H$(fo\F) < oo, then 

4;H*(f(t)\F) ^ (t ^ 0) (5.2) 

where f(t) stands for the (unique) solution to (15 -lb - Moreover, still arguing as in 11281 . o«e proves 
that if moreover / (l + v 2 + | log/o(v)|) /o(u) du < oo, f/iera 

lim / |/(w,t) - F(v)\ dv = 0. 

6. Regularity of the steady state 

In this final section, our aim is to establish the existence of some smooth stationary solution 
to (11- lb - Namely, adopting the strategy of [32, Section 4.1], we prove 

Theorem 6.1 (Regularity of stationary solutions). There exists a stationary solution F to the 
Boltzmann equation 

tQ{F,F)+C(F) = 

that belongs to C°°(R 3 ). 

We shall follow the same lines of 071 Theorem 5.5] and 021 Section 3.6], from which we de- 
duce the exponential decay in time of singularities and thus the smoothness of stationary solutions. 
This proof needs the following ingredients: 

i) The stability result already proved in Proposition 14.71 

ii) An estimate on the Duhamel representation P2l Proposition 3.4] of the solution to (11.11 ) (see 
Proposition I6T21 . 

iii) A result of propagation of Sobolev norms (see Proposition I6.3I ). 

Let us first extend the regularity estimate of l32l Proposition 3.4] to our situation. For any 
feL\ let 

^(f)( v ) = T (I " I * /) ( v ) + v i v ) = T I \ v ~ w \f( w ) dw + v (v). 

It is easy to see that, for /q £ L\, the unique solution f(t) to (11.11 ) is given by the the following 
Duhamel representation: 

f(v, t) = f (v) e" E M W '«) ds + f (r Q+ (/, /) + £+ (/)) (v, s) e" £ W)M dr dfi 

( J ° (6-1) 

= f (v)G(v,0,t) + [ (rQ + (f,f) + C + (f))(v,s)G(v,s,t)ds 
Jo 

where we set 

G(v, s, t) = exp (- J S(/)(«, r) dr^j s < t, v E M 3 . 



GRANULAR GASES DRIVEN BY INELASTIC SCATTERING 



25 



Proposition 6.2. There are some positive constants Crjuh) K such that for any k G N and ij ^ 

we have 



\\f (-)G(;0,t)\\ H k + i <C Duh e- Kt \\M\ Hk+1 ( sup ||/(-,r)||^ + sup \\f(-,r)\\ k + 3 



and 



I G(-,s,t){r Q+(/,/) + £+(/))(-,*) da 
J o 



(6.2) 



(6.3) 



< (7 Duh sup ||/(-,r)|| ft + sup ||/(-,r , 

Proof. The proof is quite similar to [37 , Proposition 5.2]. Here, for simplicity we have done it for 
natural k, although it is simple to generalize it to k > by interpolation. Precisely, for any / G L 1 
define 



L (f)i v ) = \v - w\f(w)dw. 
7r3 

It is clear that 

E(f)(v)=L(rf + F 1 )(v), V/gL 1 . 
Now, according to |23l Lemma 4.3], for any given k ^ and any 5 > 3/2, the linear operator 

L : Hg — > W^ 1 ' 00 



is bounded, i.e. for any 5 > 3/2 and any k ^ 0, there exists Cfc,5 such that 

ll L (#) II w w,°° ^ ^mII^IIh^ v f G -^5- 

Let us fix now k G N and J > 3/2. Since Fj G H% due to Assumption 12.11 one deduces that 

||E(/)|| W ,+ 1 , 00 <C||/ + F 1 ||^, V/Gtf£ 



where, as in the rest of the proof, we shall denote any positive constant independent of / and 
possibly dependent on Fi by C. Setting 

F(v,a,t) = [ £(/)(«, r)dr, 

J s 

one sees that 

\\F(;s,t)\\ w k+i,oo ^CVt^(J*\\f(-,r)\\ 2 H ,dr\ ' + C(t - s) {{F^, O^s^t. 

Now, since L(g) ^ for any g ^ 0, according to Assumption 12.11 and (12.111 ). we see that there 
exists some constant \ > such that 

^ X(Fi)(u) > x, V/ G L 1 , / ^ 0, to G M 3 . 
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By taking the successive derivatives of G(v, s, t) = exp(— F(v, s, t)), one gets as in [ 37 , Proposi- 
tion 5.2] 



ft \ (fc+l)/2 

J s ||/(-,r)||^drJ +(t- S )||F 1 ||^ + l 



<:Ce~ K ^ (l + sup \\f(-,r)\\ k+ A 



(6.4) 



for some < K < x- Then, we shall use the following estimate [37, Lemma 5.3] that allows to 
exchange a time integral and a Sobolev norm: 

/ Z(-,s)ds ^-^=( [ e A( *~ s) ||Z(-,s)||^-ds) , VA>0, Vf,r£l. 
Jo hi vA \Jo 1 J 

As a consequence we have for any k ^ 0, 

f (rQ+(/, /)+£+(/)) (- )S ) G(-,s,t)ds 
Jo 



^ C 



I e K ^ (\\rQ + (f,f)(;s)\\ 2 KVi + \\C + (f)(; \\G(; 



ds 



1/2 



Recall now the so-called Bouchut-Desvillettes-Lu regularity result in Propositions 13 . 1 1 and 13 .21 



\Q + (f,f)\\ H ^^C 

17+1 



ft + 
V+3 



»7+3 



and 



1 "-"17+1 



'Hlfl* lU llfl* + ll F HlLi lU Hii 



?7+3 77+3 ?7+3 ^7+3. 

Arguing now as in ||37l Proposition 5.2] and using the estimate (16.41 ) with the choice 5 = rj + 3, 
we get 



[ (rQ+(/, /) + £+(/)) (-, S ) G(-,s,t)ds 
Jo 

'I 

ll/MIlk e- 2 ^)fl+ sup H/C-.r)^ 1 V d S 



1/2 



^Cmax( sup ||/(-,r)||^ , sup ||/(-,r)|fe 3 
which proves ( I6.3I ). The proof of ( 16.21 ) is similar. 



□ 



A direct consequence of the previous result together with the uniform I? bounds is the uniform 
in time propagation of Sobolev norms. The proof is carried on exactly as in [32, Proposition 3.5]. 
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Proposition 6.3. Let Fi satisfy Assumption \2.1\ Let fo 6 L\, fo ^ with unit mass and let f 
be the unique solution of the Boltzmann equation dl.ll ) in C(R + ; Ljj,) D Lg) associated 

with /q. Then, for all s > a«<i 77 ^ 1, f/zere ex/sto > smc/z ?/ia/ 

f £H s +w sup||/(.,t)||« <+oo. 

The previous ingredients allow to proof the following theorem, see 11371 Theorem 5.5] for the 
proof. 

Theorem 6.4 (Exponential decay of singularities). Let fo € L\ n L 2 w/f/j wmY mass and let f 
be the unique solution of the Boltzmann equation (11.11 ) in C(R + ; L\) n Lg) associated 

with Jq. Let Fi satisfy Assumption \2.1\ Let s ^ 0, q ^ be arbitrarily large. Then f splits into 
the sum of a regular and a singular part f = fn + fs where 

( SUp \\fR(t)\\ H s nL l < +OO , f R > 

\ 3A>0 : \\fs(t)\\Ll=0{e- Xt ). 

Proof. The proof is easily adapted from ll32l Theorem 3.6] since the ^-stability result (Proposi- 
tion !4.7b . the Duhamel representation (Proposition l6.2l ). the uniform propagation of Sobolev norms 
(Proposition 16.31 ) allow to adapt directly 1371 Theorem 5.5]. □ 

Finally, Theorem 16.41 allows to prove the main Theorem 16. II 
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